INTERVAL METHODS FOR 
FMR SPECTRA SIMULATION 

Marek W. Gutowski 
gutow@ifpan.edu.pl 



Part of author's statutory activity in: 
Institute of Physics 
Polish Academy of Sciences 
Warsaw 



Marek W. Gutowski, M^dralin, August 29th - September 1st, 2002 



What is interval calculus? 



A branch of (numerical) mathematics, which makes 
possible to evaluate ranges (bounds) of algebraic 
expressions over finite domains, not just for specific 
values of their parameters. 

The results obtained on different computers may differ, 
i.e. the bounds are usually overestimated, but always 
include the true result. 



Examples: 

([-2,3]f = [0,9] 



-2, 3] + [0, 1] = [-2, 4] 



but, surprisingly 



[0, l]-[0, 1] = [-!,+!] 
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What is FMR? 



Ferromagnetic resonance is the resonance absorption 
of electromagnetic radiation by magnetic bodies placed 
in an external magnetic field. 

The magnetization vector of the ferromagnetic body, 
which is located in an external magnetic field, is freely 
precessing around the field vector with some frequency. 
The precession is damped and spontaneously decays. 
The relation between the field strenght and the 
precessing frequency is simple: 

UJ = 7He// 

where Hg// is an effective field acting on the body, 
and 7 is a constant known as gyromagnetic ratio. 

He// includes: external magnetic field, 
demagnetizing field, due to the sample's geometry 
(shape), and anisotropy field. 
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Two types of measurements 



There are two possible kinds of experiments: 

1. with fixed external magnetic field and variable 
frequency, and 

2. with fixed (microwave) frequency and variable 
external field 

We will talk about the spectra obtained with the 
second method. Here the resonance absorption is 
observed only for some magnitudes of an external 
field, called resonance fields. Due to the presence of 
various internal fields, the observed resonance fields 
will be generally different for different orientations of 
the sample in the external field. 
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The resonance condition 



The resonance condition is usually written in the form 
of the well known relation: 




d^E d^E ( d^E 



80'^ dip^ \d0 dip 



(1) 



where the second derivatives of the free energy 
^(He^t5 M, . . .) are taken at equilibrium position 
of the magnetization vector M, and (O^cp) describe 
the orientation of this vector in the polar reference 
frame. 

In order to find the equilibrium position of M, one 
needs first to solve the system of equations 

dE/de = 
dE/dip = 

and making sure that its solution(s) are indeed at the 
minimum of the free energy, i.e. that the r.h.s. of (1) 
is strictly positive. 
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Outline of classical calculation method 



1. fix the orientation and magnitude of the external 
field 

2. find numerically the equilibrium position of 
magnetization vector M and verify that it is stable 

3. calculate the resonance frequency uj 

4. if a; coincides with the frequency used in experiment 
then we have found the resonance field, otherwise 
the calculations should be repeated for other value 
of H^xt (with the same orienatation). 

The equilibrium position of M may be hard to find 
(and thus inexact), since even in amorphous, i.e. non- 
crystalline, samples, usually M and Hg^t ai'e NOT 
parallel. 
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Interval method 



The list of 3D boxes x A6^ x AH^xt is systematically 
reviewed, starting from the single initial box 

[0, 27r] X [0, TT 

For each box in succession the series of tests are 
applied, leading either to elimination of the box from 
the list or to its splitting into two (smaller) offspring 
boxes. Failing any of the tests below (answer: NO) 
eliminates the box from list 

• does the interval dE/dip contains zero? 

• does the interval dE/dO contains zero? 

• does the interval describing r.h.s of relation (1) 
contains positive numbers? 

• does the interval uo, computed from relation (1) 
contains the experimental frequency ooexp^ 



X [0, HjYiax\ 
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Interval method — some details 



Boxes not failing applied tests remain on the list. The 
largest of them is then selected and divided into two 
parts, each of which are tried again. We continue this 
procedure until the list is empty or contains only small 
boxes, i.e. in our case Ayp = A9 < 2 • 10~^ rad ^ 
10~^ degree, AH = 0.005 Oe. The maximum length 
of list is usually close to 200, but occasionally, for 
some "difficult" orientations, it exceeds 2000. On exit, 
the small neighboring boxes are 'glued' together, if 
necessary. In rare cases, for some directions of H, this 
procedure leads to higher inaccuracy in determining 



res ■ 



For amorphous wire, depending on orientation of Hext, 
zero, one, two or even more boxes (resonance fields) 
are returned, see figures. 

Typical running time, for Oext = 0°, 2°, 4°, . . . 180° is 
around 12 min. on a 100 MHz PC. 
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Advantages of interval method 



• no resonance field is ever missed 

• connplete elinnination of nunnerical inaccuracies 

• the equilibrium positions are calculated exactly, 
without any simplifications, even when the 
anisotropies are quite complicated 

• reliable replacement for other methods 

• the method may be easily extended to reliably 
reconstruct the values of unknown anisotropy 
constants and other material parameters, together 
with their uncertainties, from experimental data, 
thus replacing the usual trial-and-error procedures. 
In this case no classical counterpart — other 
than guessing — exists. 
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Disadvantages of interval method 



• neglible or non-existent knowledge of interval 
methods among the practitioners in the field and, 
generally, among physicists at large^ 

• increased requirements for the raw computing power 
(moderate) 



To find more on interval methods click on the URL: 

http: //www. cs .utep. edu/interval-comp/ 
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List of figures 



All figures simulated with: uJexp = Stt x 9.243 GHz, 
g = 2.00 and AirMg = 6400 Gs. Anisotropy constants 
Ku and are given in lO^xerg/cnn'^, angles (x-axis) 
in degrees, and resonance fields (y-axis) - in Gs. 

1. No anisotropy at all, = = 



2. 


Ku = 


-1.11, 


K^ = 





3. 


Ku = 


+0.20, 


= 


-Ku 


4. 


Ku = 


+2.00, 


K^ = 


+4.41 


5. 


Ku = 


+ 1.59, 


K4 = 


+2.85 


6. 


Ku = 


+ 1.19, 


K4 = 


-2.38 


7. 


Ku = 


-0.19, 


K4 = 


+2.90 


8. 


Ku = 


-1.19, 


K4 = 


+7.16 


9. 


Ku = 


-1.39, 


K4 = 


+3.18 


10. 


Ku = 


-5.96, 


K4 = 


-1.19 


11. 


Ku = 


-7.16, 


K4 = 


+4.77 


12. 


Ku = 


-1.39, 


K4 = 


+2.38 
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